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ABSTRACT 

We consider a differentially rotating, 2D stellar disk perturbed by two steady state spiral density 
waves moving at different patterns speeds. Our investigation is based on direct numerical integration 
of initially circular test-particle orbits. We examine a range of spiral strengths and spiral speeds 
and show that stars in this time dependent gravitational field can be heated (their random motions 
increased) . This is particularly noticeable in the simultaneous propagation of a 2-armed spiral density 
wave near the corotation resonance (CR), and a weak 4-armed one near the inner and outer 4:1 
Lindblad resonances. In simulations with 2 spiral waves moving at different pattern speeds we find: 
(1) the variance of the radial velocity, <jr 2 , exceeds the sum of the variances measured from simulations 
with each individual pattern; (2) a R can grow with time throughout the entire simulation; (3) a R is 
increased over a wider range of radii compared to that seen with one spiral pattern; (4) particles diffuse 
radially in real space whereas they don't when only one spiral density wave is present. Near the CR 
with the stronger, 2-armed pattern, test particles are observed to migrate radially. These effects take 
place at or near resonances of both spirals so we interpret them as the result of stochastic motions. 
This provides a possible new mechanism for increasing the stellar velocity dispersion in galactic disks. 
If multiple spiral patterns are present in the Galaxy we predict that there should be large variations 
in the stellar velocity dispersion as a function of radius. 
Subject headings: stellar dynamics 



1. INTRODUCTION 

The observed correlation between the ag es and veloc- 
ity dispersions of solar neighborh ood stars ijEggen et al.l 
119621 lDennislH966l lWielenlll977tl has been a sub j ect of 
study since the work of iSpitzer fc Schwarzschildl l|1951I 
1953). They established that scattering of stars from 
initially circular orbits into more eccentric and inclined 
orbits was a likely explanation for the increase in velocity 
dispersion, er, with age t. They suggested that massive 
gas clouds (then undetected) were responsible. Molecu- 
lar cl ouds were thought to be the sole scatterin g agents 
(e.g.. lMihalas fc Binnevll981fl until lLacevI l|1984ft showed 
that the observed ratio of the dispersion in the direction 
perpendicular to the Galactic plane and that toward the 
Galactic center, <j z /(Tr, was too low to be consistent with 
the predictions from this scattering process. There was 
also a discrepancy between the predicted and then mea- 
sured relation between a and t: if a oc t Q , models pre- 
dicted a ~ 0-20j whil e observations of the time yielded 
a ~ 0.5 (| WielerJI 1 977|) . This resulted in the development 
of models that incorporated the h eating of the stellar 
disk from transient spiral structur e ([Barbimis_fe_Woltjer| 
19671 ISellwood fc Carlberel Il98i ICarlberg fc Sellwoodl 
1985) in addition to scattering from molecular clouds 
ll Jenkins fc Binnevlll990llJenkinslll992|h Other proposed 
models for the heating of stars includ e scattering by 
halo black holes llLacev fc Ostrikerll985(l or dark clusters 
l|Carr fc Lacevl 119871). giant molecula r cloud s and halo 
black holes lIHanninen k. Flvnnll20ql I2004D . and infall 
of satellite galaxies ijVelazauez et al.lll999j) . 

We define a steady state spiral pattern to be a spi- 
ral perturbation moving at a constant angular velocity 
with a fixed amplitude and a fixed pitch angle. An 
individual steady state spiral arm pattern is not ex- 
pected to increase the velocity dispersion of a stellar 



population (e.g.. iBinnev fc Tremaind|1987D. though this 
may be violated at resonances (e.g., lContopoulosllT981t 
ICombes et all 119901 |Fu3 12001J) . Consequently, models 
used to explain the age/dispersion relation rely on tran- 
sient spiral structure. In the solar neighborhood, at a 
distance of 8 kpc from the Galactic center and with an 
angular rotation speed of 220 km/s one rotation period 
for a star corresponds to about 0.25 Gyrs. Thus, only 
about 40 rotation periods have taken place during the 
lifetime of the Milky Way disk. Because of the small 
number of periods, it is difficult to apply a diffusive the- 
o ry to the proble m of disk heating. 

Quillen (2003) pointed out that when perturbations 
at two pattern speeds, such as a bar and a spiral den- 
sity wave are present in the disk, the stellar dynamics 
can be stochastic, particularly near resonances associ- 
ated with one of the patterns. While this paper focused 
on spiral structure near the Galactic bar's outer Lind- 
blad resonance, here we consider the more general possi- 
bility that there could be more than one wave present 
in the disk. There is also observational evidence for 
this in other galaxi es, such as asymm etries in the spi- 
ral structure (e.g., iHenrv et al.l l2003|) . By expanding, 
galaxy images in Fourier components, lElmegreen et alJ 
( 1992) noted that many galaxies exhibit hidden three 
armed components and suggested that multiple spiral 
densit y waves can propa gate simultaneously in galaxy 
disks. iRix fc Riekel ljl993f) . in their Fourier expansion of 
the near-infrared images of M51 saw both m = 1 and 
m = 3 components as well as the dominant 2 - armed 
structure. In their recent work iNaoz &: Shavivl (2005) 
found evidence for multiple spiral patterns in our home 
galaxy By studying the nearby spiral arms they find that 
the Sagittarius-Carina arm is a superposition of two fea- 
tures, moving at different pattern speeds. Clumps con- 
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taining ol d stars in the Solar neighborhood velocity dis- 
tribution iDelmenl OjlQM) ■ may be interpreted in terms of 
multiple stellar densi ty waves traveling in ou r vicinity in 
the Milky Way disk (pe Simone et alJl2004 - Theoreti- 
cally this situation corresponds to scenarios tha t include 
multi ple and tran sient spiral density waves fe.g.. lToomrel 
ll98lUFuchsi200H) rather than thos e that focus on a single 
dominant quasi-steady mode fe.g.. lLowe et alJll994|) . 

In this paper we are interested in the effect on the ra- 
dial velocity dispersion and on the orbits of stars caused 
by the simultaneous propagation of two steady state 
spiral patterns in a combination of (strong) 2-armed 
+ (weak) 4-armed, where the second pattern moves at 
an angular speed different than the first one . Thi s 
situation has been motiv a ted bv | Le pine et ail l|2001|) . 

200l . INaoz fc Shavivl <|2005|) . "and 
where a combination of 2- 
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armed and a 4-armed spiral waves has been considered 
in the Solar Neighborhood. 

We study the kinematic and spatial evolution of stars 
in a 2D disk by following the trajectories of non-self- 
gravitating particles under the gravitational potential of 
the background (disk and halo) and the 2 spiral den- 
sity waves just described. We explore a range of pattern 
speed combinations and look for an increase in the radial 
velocity dispersion of particles. 

Our simulation model is described in fJ21 We describe 
our simulations in <J3"3S] A summary and discussion fol- 
lows in SjHl 

2. NUMERICAL MODEL 

2.0.1. Notation and units used 

For simplicity we work in units in which the galactocen- 
tric distance to the annulus in which stars are distributed 
is i?o = 1, the circular velocity at Rq is Vq = 1, and the 
angular velocity of stars is flo = Vo/Ro — 1- Actually, 
since we assume a flat rotation curve throughout the pa- 
per the initial circular velocity is Vq = 1 everywhere. 
One orbital period is 2ttRq/Hq — 2n. In our code time is 
in units of 1 /fi whereas in the figures it is given in units 
of orbital periods at Rq . The velocity vector of a star is 
(u, v), where u, v are the radial and tangential velocities 
in a reference frame rotating with Vq . Consequently, the 
tangential velocity of a star in an inertial reference frame 
isVo + v. 

In the definitions that follow the subscript a refers to 
the primary, 2-armed spiral wave; similarly the subscript 
b refers to the secondary, 4-armed one. The azimuthal 
wavenumbers of each spiral wave, m a , rrib, are integers 
corresponding to the number of arms. The amplitudes of 
the spiral wave gravitational potential perturbations are 
denoted as e a , £f>- The spiral pattern frequencies of each 
spiral wave are £l s ,a and £l s ,b- The parameters a a , ctb are 
related to the pitch angles of the spiral waves, p a , pb, as 
a = mcot(p), negative for trailing spirals with rotation 
counterclockwise. In this paper we only consider trailing 
spiral waves. 

We write the Hamiltonian of a star as H = Ho + Hi 
where H is the unperturbed, axisymmetric part and Hi 
is a small perturbation due to the 2-armed, 4-armed, or 
both spiral waves. We write Hi = $ a + where $ a , $b 
are the gravitational potential perturbations due to each 
pattern (see next section). 

In general we refer to the radial velocity dispersion as 



o~r. However, we also define 3 other symbols referring to 
the type of perturber we use in the Hamiltonian. Con- 
sidering the above described form of Hi , if <I> a = then 
(tr — > OR.b (only 4-armed pattern present); if $b = 0, 
or — > o-R,a (only 2-armed pattern present); finally, if 
both 3>a,3>b 7^ then or — > o-R,ab (both spiral waves 
present). 

2.1. Equations of motion 

We consider the 2D motion of a test particle in the 
mid-plane of a galaxy in an inertial reference frame. In 
plane polar coordinates, R, 4>, the Hamiltonian of a star 
can be written as 

H{R,(j) : p R ,p^t) = H {R,p R ,p 4> ) + Hi{R,<j),t) (1) 

The first term on the right hand side is the unperturbed 
axisymmetric Hamiltonian 

= A , A. 

2 2R 2 



H (R,pr,'j 



MR) 



(2) 



where we assume the axisymmetric background potential 
due to the disk and halo has the form $o(R) = Vq \og(R), 
corresponding to a flat rotation curve. Energy and an- 
gular momentum for Hq are conserved. 

In the case of one periodic perturbation there is still a 
conserved quantity in the reference frame rotating with 
the pattern. That is the Jacobi integral, J = E — L£l s , 
where E is the energy of the particle, L is its angular 
momentum, and f2 s is the pattern angular velocity. How- 
ever, upon adding a second spiral density perturbation 
H does not admit any isolating integrals. 

The perturbation is Hi(R,(/),t) = <5> a (R,<p,t) + 
$b(R,4>,t) where the 2 terms on the right hand side are 
the gravitational potentials due to each spiral wave. We 
define the subscripts a, b to denote the 2- and 4-armed 
spiral patterns, respectively. The perturbations are dis- 
cussed in the next subsection. 

Hamilton's equations applied to Eq. ^give us the equa- 
tions of motion 



(3) 
(4) 



In the first of these equations the primes denote partial 
derivatives with respect to R. 

2.1.1. Perturbation from spiral structure 

We treat the spiral patterns as small perturbations to 
the axisymmetric model of the galaxy by viewing them as 
quasi-steady de nsity waves in a ccordance with the Lin- 
Shu hypothesis ijLin et al.lll969|) . Each spiral wave grav- 
itational potential perturbation is expanded in Fourier 
components as 

$(R, <f>, t) = e m exp [i(a InR- m(<p - fl s t - <j) m ))] 

m 

and its corresponding surface density is 
E(i?, cj>, t) = S m exp [i(a In R - m(<j> - Q s t - m ))] 

(6) 
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We assume that the amplitudes, e m ,£ m , and the pitch 
angles are nearly constant with radius. The strongest 
term for a 2-armed spiral is the m — 2 term and similarly 
m = 4 for a 4-armed structure. Thus only the terms 
corresponding to those m values are retained. In this 
paper we only consider 2- and 4-armed spiral waves, thus 
in our notation m a — 2 and m& = 4. Upon taking the 
real parts of Eqs. |5] and H3 the perturbation due to the 
2-armed spiral density wave becomes 



$ a (R,<p,t) = e a cos (a a In R - m a ((f)- 



<t>0,a))- 



Changing the subscript a to b in the above equation gives 
the 4-armed perturbation. The direction of rotation is 
with increasing and a a , a% < ensure that each pattern 
is trailing. 

The recent study of iValled (|2p05) provides a good 
summary of the many studies which have used obser- 
vations to map the Milky Way disk. Cepheid, HI, 
CO and far-infrared observations suggest that the Milky 
Way dis k contains a 4-armed tig htly wound structure, 
whereas [Drimmcl fc Soersreil 1)200 1|) have shown that the 
near-infrared observations are consistent with a domi- 
nant 2-armed structure. The dominant 2-armed and 
weaker 4-armed structu re was previously proposed by 
lAmaral fc Leninel l|1997f l. We adopt the same configu- 
ration in our model - a primary 2-armed and a weaker 
4-armed spiral wave perturbations. What is different 
and new in our model is the assumption that there is 
a non-zero relative angular velocity between the 2 spiral 
patterns. This introduces an additional parameter - the 
pattern speed of the secondary spiral wave. 

In general, an individual stellar orbit is affected by the 
entire galaxy. However, if tight winding of spiral arms is 
assu med only local gr avitational forces need be consid- 
ered. iMa et all l|200ff) find kR for several spiral galax- 
ies of various Hubble types to be > 6, thus the tight- 
winding, or WKB approximation, is often appropriate. 
In the above expression k is the wave- vector and R is the 
radial distance from the Galactic Center (GC), related to 
the pitch angle p through cot(p) = \kR/m\. This gives 
a in terms of k as a f» \kR\. In the WKB approxima- 
tion the amplitude of the potential perturbation Fourier 
component is related to the density perturbations in the 
following way 

—2i:GYiQSaRo . , 



ijBinnev fc Tremaind 11987^) . The above equation is in 
units of Vq . E a is the amplitude of the primary, 2-armed- 
spiral-wave mass surface density and S a » £„/£()• A 
similar relation holds for the second spiral pattern. 

2.1.2. Resonances 

Of particular interest to us are the values of fl s which 
place the spiral waves near resonances. The Corotation 
Resonances (CR) occurs when the angular rotation rate 
of stars equals that of the spiral pattern. Lindblad Res- 
onances (LRs) occur when the frequency at which a star 
feels the force due to a spiral arm coincides with the 
star's epicyclic frequency, k. As one moves inward or 
outward from the corotation circle the relative frequency 
at which a star encounters a spiral arm increases. There 



are 2 values of R for which this frequency is the same 
as the epicyclic frequency and this is where the Outer 
Lindblad Resonance (OLR) and the Inner Lindblad Res- 
onance (ILR) are located. 

Quantitatively, LRs occur when fl s = Slo i njm. The 
negative sign corresponds to the ILR and the positive - 
to the OLR. Specifically, assuming a flat rotation curve, 
for the 4:1 ILR Vl s = 0.65fi , for the 4:1 OLR n s = 
1.35O 0) for the 2:1 ILR tt s = 0.3Q Q , and for the 2:1 OLR 
Q s = 1.7Qq. Corotation of each spiral pattern occurs at 

= "0- 

2.2. Simulation method 

To investigate the effect of a second spiral density wave 
on the orbits of stars, we ran simulations of test particles 
in a disk with two spi ral density waves. As done by 
iDe Simone et alJ <)2004[) the amplitude of each wave was 
varied, however the pitch angle and pattern speed were 
held fixed during the integrations. We are restricting this 
initial study to this type of spiral arm growth (no swing 
amplifier mechanism). 

2.2.1. Initial conditions 

Positions of test particles (stars) were chosen randomly 
with a uniform density distribution in both radial and 
azimuthal directions. Thus, we distribute stars in an 
annulus of an inner and outer radii Rq — AR, i?o + AR, 
where for most of the runs Ai? is 0.3. All stars were given 
the same initial velocity (u, v) = (0, 0) consistent with a 
flat rotation curve and an initially cold stellar disk. 

We now describe the growth of the spiral density waves. 
The amplitudes, e a ,e{, are zero before to, grow linearly 
with time at < t < t\ and stay constant at e a — e maXja , 
£fc = tmax,b after t = t\ = 4 rotation periods. 

Simulations were run for 40 rotation periods, which 
corresponds to about 10 Gyr. Heating rates were es- 
timated calculating the radial velocity dispersion (the 
standard deviation of the radial velocity component u), 
every rotational period. This was done for all stars in 
the annulus (R Q - Ai?,i? + AR) where Ai? = 0.05. 

The standard deviation was computed as the square 
root of the bias-corrected variance: 



i 



(9) 



where n is the number of stars in the bin (the annulus 
of width Ai? = 0.05 centered on i? , described above) 
and u = Q^iLi u i)l n - Table 1 shows the parameters for 
which simulations were run. 

2.2.2. Numerical accuracy 

Calculations were performed in double precision. We 
checked how energy was conserved in an unperturbed test 
run. The initial energy was compared to that calculated 
after 40 periods (« 10 Gyr) and the relative error was 
found to be \ AE/E(G)\ < 8.48 x 10~ 14 . We also checked 
the conservation of the Jacobi integral J in the presence 
of one periodic perturbation. The relative error found in 
this case was |AJ/J(0)| < 3.11 x 10~ 12 . 

3. ONE SPIRAL PATTERN 
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To tell the difference between the effect of one spiral 
pattern and both patterns, we first ran comparison 
simulations with only one spiral pattern present. In 
each run we integrated 6000 particles. Parameters used 
for integrations are listed in Table 

In FigQ]we plot the radial velocity dispersion, OR. a as 
a function of time for the primary 2-armed spiral wave 
only. In all panels the amplitude of the gravitational 
perturbation is kept fixed at e a — —0.006 once the wave 
is grown, and the pitch angle function is a a — — 6 which 
corresponds to p a ~ —18°. Each panel in Fig. Q shows 
the increase of the radial velocity dispersion with time at 
a different pattern speed, Q s , a - 

Logarithmic spirals are known to be self-similar. Due 
to this fact, changing the angular frequencies of the spi- 
ral waves is equivalent to changing the radius. Since 
we assumed a flat rotation curve, Vq = 1 everywhere 
and the angular velocity of stars varies with radius as 
Q = Vq/R~ 1/R. On the other hand, since the pattern 
speed, Q s does not vary with R, i.e., the pitch angle is 
constant, a spiral rotates as a rigid body. It follows that 
in the units used in this paper we can just substitute R 
for fl St a in Fig. ^ placing Rq at £l Sya = 1.0. We note 
that since stars in the inner parts of a flat-rotation-curve 
galaxy complete one rotation period faster than stars in 
the outer parts, 40 rotation periods at i?o translates as 
a larger number of revolutions in the inner parts and a 
smaller number in the outer parts. Assuming a corota- 
tion radius near Rq and a galactic radius of ~ 2Rq, Fig. 
Q gives the heating as a function of radius in the galaxy 
for one pattern. The top left hand panel shows pattern 
speeds placing Rq near the 2:1 ILR. The lower right hand 
panel shows pattern speeds placing R near the 2:1 OLR. 

As the spiral amplitude, e a , reaches its maximum value 
the radial velocity dispersion levels off. As predicted 
by many previous theoretical investigations, there is no 
heating after the wave is grown. As expected, large 
radial velocity dispersion values are found only near 
the 2:1 Lindblad resonances (0,<,. n , = 0.3 , 1.7). This 
is consistent with iLvnden-Bell fc Kalnaisl l)1972t ) who 
showed that steady spirals heat stars only at the Lind- 
blad resonances. Note that at and around corotation 
(fis.a = 1.0) random motions of stars were not increased. 

The effect of a weak 4-armed spiral pattern on the 
radial velocity dispersion of stars is shown in Fig. |2 
for 3 values of the gravitational potential perturbation 
amplitude, Different line styles show tb — —0.001 
(dotted), e b = -0.002 (dashed), and e b = -0.003 (solid). 
As in the case of the 2-armed spiral wave, relatively 
large ctR.b values are observed only near the first order, 
4:1 LRs (fi s ,6 = 0.65,1.35). Substituting R for Q s>b in 
Fig. Ogives <jr=<jr(R) in the range (R = 0.5, R = 1.6). 

In both the 2- and 4-armed spiral structures the ini- 
tial rapid increase ends when the spiral amplitudes attain 
their maximum values (t « 4 rotations) and are fixed for 
the rest of the runs. This observation will later be con- 
trasted with the effect of the simultaneous propagation 
of 2 spiral waves (see H4.2|) . 

4. TWO SPIRAL PATTERNS GROWN SIMULTANEOUSLY 



We next ran simulations with two spiral density waves 
moving at different pattern speeds. The spiral ampli- 
tudes were grown simultaneously over 4 rotation periods 
as in the case of only one spiral wave. 

4.1. Heating in (T2 Sj0; Q, s ^)-space 

First we carried out a set of integrations with 1000 
particles going through all the possible combinations of 
the two pattern speeds in the range (0.1,2.0) for the fi Sja 
and (1.5,1.6) for f2 s ^, in steps of 0.1, and in units of 
the angular rotation rate, Slo at Rq. The chosen ranges 
of £l Sl a and f2 Si b ensure that we cover the Lindblad and 
corotation resonances of each spiral wave. The results are 
shown as the (ft s ,a, f2 S! b)-space contour plots in Figures 
[3]and0J described below. 

Horizontal axis is Q s , a and vertical axis is Si f, for each 
panel. Each panel corresponds to a different secondary 
spiral wave amplitude, e b = -0.001, -0.002, -0.003. The 
primary pattern amplitude is constant for all panels at 
e a = 0.6. Spiral amplitudes were grown simultaneously 
in 4 rotational periods and were kept fixed after that. 

In Fig. 13 contours show 

ACTR = (T R:ab - \JoR,a 2 +<JR,b 2 - (10) 

Here OR. a and <jrj, are the radial velocity dispersions as 
a result of only the primary or only the secondary spiral 
patterns, respectively. These were time averaged over 40 
galactic rotations, excluding the time during which wave 
amplitudes grew (initial 4 periods). Similarly, OR^b is 
the time average of a simulation run with both patterns 
perturbing the test particles. Consequently, A/jr gives 
the additional heating caused by the combined effect of 
the two perturbers as opposed to the sum of the heat- 
ing that would have been produced if each spiral wave 
propagated alone. 

The dashed lines show the locations of corotation and 
Lindblad resonances due to each of the two spiral wave 
perturbations. Next to each dashed line we labeled the 
corresponding resonances by denoting the corotation res- 
onance with " CR" , and the outer and inner Lindblad 
resonances with "OLR" and "ILR". 

Grayscale for A/jr was normalized to the maximum 
value of the three panels allowing a proper comparison 
among the panels. Darker regions corresponds to larger 
values of A<jr. 

We see that Agr > for most combinations of an- 
gular velocities which implies that spiral waves interfere 
and produce additional heating. In all the panels maxi- 
mum values of Acfr are attained mainly along the lines 
ft s .a—l-0 and f2 S) {,=1.0 corresponding to the CR of each 
spiral wave. A notable exception to this trend is the 
intersection of these lines where patterns do not have 
relative motion with respect to each other. As the abso- 
lute value of the secondary spiral wave amplitude, |e&|, is 
increased from 0.001 to 0.003 structure in the fl Sta ,fl St b- 
plane seems not to change significantly, however A<jr is 
increased. 

It is interesting to note the lack of heating due to the 
effect of the coupling between the two spiral waves at the 
2:1 ILR and on the outside of the 2:1 OLR of the 2-armed 
spiral wave, and at the vertical line f2 a « 0.7 which is just 
outside the 4:1 ILR with the same perturbation. 

We mentioned earlier that changing Q s and keeping 
R constant is equivalent to changing R and keeping f2 s 
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constant. Thus, for a given set of pattern speeds, a line 
at 45° slope in Fig. El predicts the variation of Actr with 
radius, R, where the intercept of the line is set by the 
relative angular frequency, AQ S = fl Sta — Q s ,b- Depend- 
ing on the given AQ S , a galaxy heated by two spiral 
density waves would exhibit different variation of the 
radial velocity dispersion with radius. For example, if 
Af2 w then the galaxy would not be heated efficiently. 
This is expected since, as mentioned earlier, steady 
spirals heat stars only at the LRs. If, on the other 
hand, Af2 rs 0.3 then the currently examined heating 
mechanism would prove quite efficient and would effect 
stars over a wide region. It is evident from the figure 
that heating takes place mainly near the CR and OLR 
of primary, and near the CR of the secondary pattern. 

Similarly to the Actr contours shown in Fig. EJwe plot 
contours of the radial velocity dispersion as the result 
of the simultaneous propagation of the two spiral waves, 
a Rt ab, in Fig. H 

Here maximum values are attained at the 2:1 ILR of 
the primary, e a = —0.006 spiral wave. We believe that in 
reality the 2:1 ILR is too strong to maintain spiral struc- 
ture for strong spirals so we regard the results observed 
in this location on the f2 Sja , f^^-plane as unphysical. It 
is possible that in reality spiral amplitude decreases as 
this resonance is reached or that spiral arms end there. 
Aside from the 2:1 ILR of the 2-armed spiral wave in 
all the panels high values of a^ab are observed around 
the 2:1 OLR of the primary spiral wave or near reso- 
nances with the primary and the secondary pattern tak- 
ing place at the same time (where dashed lines cross) . As 
in Fig. |21 an exception to the last statement is the point 
(0 SiO = 1.0, £l S: b — 1-0) for which the lowest values in 
CR,ab are reached. Again, this corresponds to both pat- 
terns placing stars at the CR and thus no relative motion 
between them. 

In the following section we chose interesting points cor- 
responding to specific combination of the pattern angular 
velocities in Fig. Eland plotted the radial velocity disper- 
sion, (jR^ab, as function of time, t, to compare to the effect 
of the single spiral waves described in Parameters for 
all runs can be found in Tabled 

4.2. Time evolution of the radial velocity dispersion at 
the CR and OLR of primary pattern 

To better understand the nature of the process giving 
rise to the structure seen in Figures E] and 01 we plot- 
ted the time evolution of or^ and */ or^ 2 + o~R,b 2 (as 
defined in the previous subsection, see Eq. I10| l. These 
dispersions are shown in Figures l5l7l as functions of time. 
We chose to consider stars near the corotation resonance 
with the primary pattern and with different values for 
the secondary pattern speed. For reference, the time av- 
erage values of Ao~r and <TR, a b are given by the vertical 
line fl s .a — 1.0 in Fig. El and the same line in Fig. 01 
respectively. 

The solid lines in Figures 15171 depict the effect of the 
simultaneous propagation of a 2- and a 4-armed spi- 
ral patterns on the radial velocity dispersion, erR, a b, as 
a function of time, t, in units of galactic rotation pe- 
riods. Figures EJ El an d differ only in the value 
of the amplitude of the secondary spiral wave, to wit, 
e b = -0.001,-0.002,-0.003, respectively. The value of 



the primary spiral wave amplitude is fixed for all 3 Fig- 
ures with e a = —0.006. The 2-armed spiral structure has 
the same parameters as that shown in Fig0 for which 
fl s .a = 1.0 (CR), whereas the 4-armed spiral wave has 
different angular velocity in each panel of Figures 15171 
The dispersion due to the latter perturbation alone, ctr^, 
was shown in Fig. El For comparison, in addition to cr^^ 
(solid line) we also plotted ^/o~R. a 2 + OR.b 2 (dashed line) 
(see Eq. EJ. 

The 4:1 resonances of the secondary spiral wave are 
seen in the lower left 2 panels (4:1 OLR) and the up- 
per middle 2 panels (4:1 ILR). Note the big dynamical 
changes at those locations for the case of a-R, a V- 

(1) OR.ab increases up to a factor of « 3 compared to 
the dispersion expected from the noninteracting spiral 
waves, \JorJ 1 + o-rj, 2 (i.e., panels with fit = 1.2, 1.5, 
Figures IST7J) : 

(2) o-R } ab grows with time. In Fig. EJat £l Sy i, = 0.6, 0.7, 1.4 
°~R,ab increases until t w 10; for f2 S) & = 0.8,1.2,1.5,1.6 
o~R,ab increases at even later times of 20-30 periods. This 
does not happen when only one pattern is present. 

(3) resonances have effect over a larger range of angular 
velocities, compared to U o-R. a 2 + UR.b 2 , and thus over a 
larger range of radii. 

As noted earlier, in the case of only one spiral density 
wave large effects are primarily observed near the 2:1 
Lindblad resonances for the 2-armed, and near the 4:1 
Lindblad resonances for the 4-armed spiral patterns dur- 
ing the growth of the perturbations only. In contrast, 
when two spiral waves at a particular combination of 
pattern speeds are present, o~R, a b keeps increasing even 
after the maximum values of the potential amplitudes, 
e a and are reached. We expect the most likely ex- 
planation for this is that the particular combination of 
pattern speeds, namely, a 2-armed spiral at the CR and a 
4-armed one near the ILR or OLR, creates stochastic mo- 
tion due to resonance overlap. This is analogou s to th e 
bar - spiral wave overlap described bv iQuillenl ( 2003), 
found to produce a stochastic region in the vicinity of 
the Galactic bar's OLR. In these conditions stars are ex- 
pected to diffuse radially in real space (see i]5.2f> . 

The dispersion ctr ab grows with time up to a specific 
value of t (saturation time) after which it does not in- 
crease anymore. It is interesting to note the decrease 
in the saturation time as the strength of the secondary 
spiral increases in Figures EE 

If the observed increase of <7R ia b with t is due to the 
presence of both the CR of the primary, 2-armed and 
the LRs of the secondary, 4-armed spiral waves, then 
one would expect that as one moves away from the 
LRs of the secondary pattern the effect on the radial 
dispersion would diminish. This is generally the trend 
observed in Figures 15171 however on both sides of the 
4:1 OLR the opposite effect is seen. At 57 Si b = 1.2,1.5 
large values of o~r a b are still evident. In fact, at the 
points (n«,„ n Slb , t)=(l,1.2, 40),(l,1.5,40) in Fig. El the 

ratio o-R,ab/ \/ &R t a 2 + o Rfi 2 attains its maximum value 
of w 3.5. These high values of Actr on each side of 
the 4:1 OLR of the 4-armed pattern were anticipated, 
considering the structure along the line Vt 8 a = 1.0 in 
Fig. El 

In Fig. El we plot the time evolution of <JR, a b and 
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\J pr.o, 2 + o~R.b 2 for Q, sa — 1.7. At this angular velocity 
of the primary pattern particles are placed at the 2:1 
OLR. The gravitational perturbation amplitude of 
the secondary wave is eb — —0.003; the primary one 
has the default value of e a = —0.006. Obviously, the 
fraction o~ R, a b / \/ o~ R,a 2 + o~R ,b 2 is much smaller than 
that observed in Figures 15171 This was expected in view 
of Figures and 0] 

It is clear from Figures I5l8l that the velocity dispersion 
in the presence of 2 patterns, o~R a b, exceeds that 
expected from each individual pattern <7R, a 2 + &R,b 2 - 
Furthermore, <JR^ a b continues to increase with time 
for some values of £l Sja , ^s,b- This contrasts with the 
predictions for heating due to transient spiral density 
waves. Additional heating occurs in the presence of 2 
waves and the interplay of the 2 waves contributes to 
additional heating after they have ceased to grow. 

We also considered the more realistic situation in which 
the second spiral density wave was grown at a later time 
than the first one. At all parameters kept the same as 
those used to produce Figures 03 H3 and0 we ran simula- 
tions in which the starting time for the secondary spiral 
wave amplitude was set to 2,4, 7, 15, and 20. We found 
that the sole effect of the delayed second spiral was to 
delay the starting point of the radial dispersion increase. 

5. ORBITS OF STARS NEAR RESONANCES IN THE 
ROTATION FRAME OF THE PATTERN 

To understand what causes the heating observed in 
Figures 15181 we will next look at the shapes of the orbits 
of test particles when subjected to only one of the 2 spiral 
perturbations and then compare those to the orbits of 
particles in the presence of both. 

The description of the particle trajectories below is 
based on movies we made for each simulation run in the 
current section, [3J Snaps from the movies are shown in 
the figures below. 

5.1. Orbits of stars subjected to one spiral wave only 

We offer a qualitative description of the orbits of stars 
near corotation with the primary spiral wave only, as seen 
in a reference frame rotating with the spiral pattern. To 
simplify this task as much as possible we integrated 1000 
particles placed at different radius R and uniformly dis- 
tributed azimuthally. Different panels in Fig. [5] show 
particles initially moving: 

(a) just inside the CR (R = 0.9), 

(b) exactly at the CR (R = 1.0), and 

(c) just outside the CR (R — 1.1). These panels are 
shown for t — 40 rotation periods for 3 separate simula- 
tions. Dashed circles show the initial distribution of the 
particles. 

Just inside of the corotation circle stars are moving 
faster than the spiral pattern (Fig^l a). In a reference 
frame rotating with the spiral, their orbits are confined, 
depending on the azimuthal angle of each star at the time 
the spiral perturbation is turned on, to one of 3 types of 
orbits. These are an oval orbit and two banana-shaped 
closed curves situated outside the initial circle, the cusps 
of which coincide with the two spiral arms. 

Just outside of the corotation radius (FigEI c ) stars 
follow similar orbits to the ones described above, except 



that because the spiral wave moves faster than the stars, 
the "banana" orbits are situated inside both the oval- 
shaped orbit and the initial circle (at Rq). Furthermore, 
the oval orbit is flattened at a 90° phase difference com- 
pared to the same type of orbit inside of the CR. 

Finally, exactly at corotation a more complicated sit- 
uation is observed (Fig|51 b). As seen in the reference 
frame of the spiral, stars move on 2 families of nested 
closed curves similar to the 2 "banana" orbits discussed 
above. Each such orbit protrudes equally radially inward 
and outward. Note, that the third, oval shaped orbit is 
not present here. Stars which happened to be at the 
arms at the time the wave was turned on, stayed there 
and the rest are not allowed to cross an arm at any time. 

We talked about the orbits of stars but what Fig. [5] 
actually shows are snapshots at a single time, t = 40 
rotations. During the entire simulation a given closed 
contour on which stars move does not change shape 
or orientation as seen in the reference frame of the 
spiral wave. Consequently, a plot showing the orbital 
trajectory of a single particle would look exactly like the 
snapshots in these pictures. However, if a contour on 
which stars move does change shape or orientation then 
it cannot be regarded as the trajectory that a single star 
would follow. 

Similarly to the CR of the 2-armed spiral wave, infor- 
mation about the orbits of stars placed at and near the 
4:1 OLR (SI;, = 1.65) of the secondary, 4-armed spiral 
wave can be obtained from Fig. ^| In all panels all 
parameters are kept the same except for the radius at 
which stars were initially distributed. The pattern am- 
plitudes are eb = —0.002, e a = and angular velocity is 
£lh = 1.65. Panels show stars initially moving on a circle 
of radius: 

R = 0.9 or just inside the 4:1 ILR (Fig. a), 

R = 1.0 or exactly at the 4:1 ILR (Fig. [TUl b), and 

R = 1.1 or just outside the 4:1 ILR (Fig. [TUl c). 

The snapshots are at time t — 40 galactic rotations. As 
we discussed at the end of the previous paragraph, the 
curve on which test particles lie at a point in time is not 
necessarily the trajectory which a single particle would 
follow in time (in a frame rotating with the pattern). It 
turns out that only Fig. ^3 (b) depicts a curve which 
changes shape and orientation as the system develops in 
time. Therefore, we do not regard Fig. ^3 (b) as the 
trajectory of a particle. 

The particle trajectories (in a frame rotating with 
the pattern) in panels (a),(c) of Fig. ^] form as soon 
as the spiral wave amplitude, eb, is grown (4 galactic 
rotations) and do not change until the end of our runs 
(40 rotation periods). However, exactly at the 4:1 ILR 
(Fig. ^3 b) the rectangular curve starts to form at 
t « 2 rotations with its vertices supporting the spiral 
structure and starts to spin clockwise in the rotation 
frame of the pattern. As it does so its vertices sharpen 
to achieve the curve seen in Fig. ^3 (b). The curve goes 
back to being oblate as the rectangle's vertices reach 
the arms again. This process is repeated with a period 
of about 9 galactic rotations. Note that star orbits 
support the spiral structure outside the 4:1 ILR fFiglTUl 
c) but are out of phase with the perturbation inside 
of it fFig llOl a). The orientation of the orbits changes 
at the resonance. It is known that strong 2-armed 
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densi ty perturbations can excite square-shaped orbits as 
well HCoutonoulos fc Orosholl ITflM IQmlleu fc Min'chTvl 
2005). However, in this case (secondary order in the 
epicyclic amplitude 4:1 ILR) the spirals are supported by 
star orbits interior to the 4:1 ILR and are not supported 
by the orbits of stars having the ir guiding rad i i exte rior 
to it. This behavior prompted iContopoulosI l|1985|) to 
propose that 2-armed spiral waves ended at the 4:1 ILR. 

To complete the picture we lastly show the behavior of 
stars when subjected to the 4:1 OLR of the secondary, 
4-armed spiral wave only. Fig. El shows a snapshot of 
the real-space particle distributions at time £=40 galac- 
tic rotations. Pattern speeds and gravitational potential 
amplitudes were the same for each panel, = 1.35 and 
e b = 0.002, respectively. The initial distribution of each 
simulation was at: 

R = 0.9 or just inside the 4:1 OLR (Fig. El a), 

R = 1.0 or exactly at the 4:1 OLR (Fig. El b), and 

R=l.l just outside the 4:1 OLR (Fig. EI c). 

As in Fig. 1101 inside and outside the resonance 
the initially circular curve did not change shape or 
orientation, thus it represents the approximately closed 
orbit of stars in a frame rotating with the pattern (Fig. 
1101 a, c). On the other hand, Fig. 1101 (b) shows a curve 
which spins as the system develops in time. 

In all the cases but Figures El (b) and El (b) , the 
initial circles assumed some particular shape at the end 
of the growth of the spiral waves and did not change 
for the remaining time. To check whether this was the 
result of an adiabatic perturbation we also tried using 
a sudden perturbation by growing the spiral in just 0.1 
period, as opposed to the default value of 4 periods, and 
found similar results. 

5.2. Orbits of stars subjected to 2 spiral density waves 
moving at different angular velocities 

In this section we attempt to look for an answer to 
the following question: What process gives rise to the 
continuous increase of the radial velocity dispersion, ctr, 
with time, t, in Figures I5l8f 

We will only look at the following combinations of pat- 
terns: 

(i) a 2-armed spiral wave at the CR (fl s ,a = 1-0) prop- 
agating simultaneously with a 4-armed one just outside 
the 4:1 ILR (fi s , ft = 0.7). 

(ii) a 2-armed spiral wave at the CR (0 SiO = 1.0) propa- 
gating simultaneously with a 4-armed one just inside the 
4:1 OLR (n Sib = 1.2). 

The sole effect of a 2-armed spiral wave, moving at a 
pattern speed which places stars at the CR, on the stars' 
orbits was shown in panel (b), Figure The trajectory 
of a star moving exterior to the 4:1 ILR (O s i, = 0.7, con- 
sidered in (i) above) is closest to the curve seen in panel 
(c), Fig. El which could be perceived as the result of a 
simulation with an initial spatial distribution on a circle 
of radius R — 1.0 and a pattern speed of Vt Si b — 0.75 
(see 3-rd paragraph of Similarly, the trajectory of a 
star moving inside the 4:1 OLR (fi s ,& = 1.2, considered 
in case (ii) above) is closest to the curve seen in panel 
(a), Fig. II II which corresponds to an initial spatial distri- 
bution on a circle of radius R = 1.0 and a pattern speed 
of fi s b = 1.25. 



Remember that in Figures T5I8I initially we distributed 
stars in the annulus (i?o — AR, R n + AR) with AR = 0.3 
and Rq = 1. Trajectories of test particles were followed 
for 40 rotations (the lifetime of the Galaxy at Rq) and 
o~r was computed every period in the annulus (Rq — AR, 
R + AR) with AR. = 0.05. Because of the differen- 
tial rotation of stars in the disk, the initial velocities of 
all test particles are fixed at {v,u) = (Vo,0) so it takes 
longer for particles at larger R to complete one period 
compared to particles at smaller R. On the other hand, 
spirals move as rigid bodies. It is obvious that particles 
moving at different radii would see the spiral structure 
moving at different pattern speed. Stars contributing to 
o~r could have come from various radii, thus responding 
to different pattern-speed spiral waves. This is a compli- 
cated process. To simplify it, considering first the panel 
for which = 0.7 in Fig. (case (i) above), we give 
the following illustration: 

(1) One can imagine the initial spatial distribution of 
particles is made of rings of stars at different radii be- 
tween i?o — AR and Rq + AR. Thus one can simplify 
the problem by looking at rings of stars. On a particular 
ring all stars feel a spiral pattern moving at the same 
pattern speed. Furthermore, a plot of the time develop- 
ment of the spatial distribution of stars would be much 
more informative than a similar plot of a simulation run 
with the set up used for Fig. panel with Q s ^ = 0.7. 

(2) Integrate 3 such rings corresponding to specific 
radii in the initial annulus distribution of Fig. El We 
chose those to be at R — 0.9, 1.0, 1.1. 

(3) Plot the spatial distribution of stars at different 
times to see whether particles, initially positioned on 
the circle of radius R, diffuse radially as was speculated 
earlier. 

Plotted in Fig. El ar e 3 steps of the 40- rotation-period 
time evolution of 3 rings of stars at the radii chosen in 
step 2 above. Rows from top to bottom show particles 
subject to 2 pattern speeds with angular velocities f2 Sja = 
1.0, fi^b = 0.7 (case (i) above) and initially positioned 
in a ring of radius: 
R = 0.9 (Fig. H2 a), 
R = 1.0 (Fig. El b), and 
R = 1.1 (Fig. El c). 

Snapshots are at times t — 10, 25, and 40 galactic rota- 
tion periods. The crosses show the position of the galac- 
tic center. 

The 3 systems in Fig. El behave very differently 
compared to the ones depicted in Figures 151 and ITUl Note 
that stars do diffuse radially. Orbits are not confined 
to narrow regions in phase space as was observed in 
Figures 1911 II We see very different dynamics in the 
trajectories of particles. In the first snap (at t=10) it 
is clear that inside the CR stars are scattered radially 
outward whereas the opposite effect is seen exterior to 
the CR. This is similar to the behavior of stars in Fig. 
El and thus must be related to the proximity of the CR. 
However, stars do not lie on a specific curve as in the 
case of only one perturbation. As one moves to panels at 
t = 25, 40 a randomization in the azimuthal direction is 
apparent. As was seen in Fig. panel with f2 s .t, = 0.7, 
saturation time was reached at around t « 10 periods. 
This is consistent with Fig. El where the real space 
particle distributions in the panels with t = 25,40 for 



each row do not appear much different. Note that, in a 
simplified way, the time development of the real-space 
distribution of stars in Fig. panel with f2 Sj b = 0.7, 
can be obtained by overlapping the 3 rows in Fig. 1121 

Similarly to Fig. El radial diffusion of stars is 
expected in the vicinity of the CR due to the 2-armed 
pattern and the 4:1 OLR of the secondary, 4-armed 
one (case (ii) above) . We show in Fig. El the time 
development of the real space particle distribution of 
3 initially circular rings of stars taken from the initial 
annulus distribution in Fig. H3 panel with Q 3; b = L2. 
The primary spiral wave angular velocity is Qb — 1.0 (as 
in Fig. IT2l placing stars at the CR, but the secondary 
pattern speed is now Oj, = 1.2, which places stars 
just inside the 4:1 OLR. We see in Fig. HO panel with 
Q s ,b — 1-2, that the extra velocity dispersion coming 
from the correlated spirals is about 3 times that due to 
the individual ones, i.e., Aer^ « 3 (see Eq. El> whereas 
in the fl s .b = 0.7 panel of the same figure Actr w 2. 
We consequently expect stronger effect on the radial 
diffusion in the case of a secondary spiral wave moving 
at f2 Si 6 = 1.2. This is indeed the case as Fig. 1131 shows. 

The effect of the addition of a second pattern at a 
different pattern speed than the first one is the observed 
radial diffusion in Figures El and El For this to 
work the 2 spirals must be correlated so that particles 
encounter each spiral wave in a systematic manner. The 
radial diffusion of stars is related to the increase of the 
radial velocity dispersion, an, with time, seen in Figures 
15181 This is most prominent around the CR or the OLR 
of the 2-armed spiral wave and the OLR and ILR of the 
secondary, 4-armed one. 

6. SUMMARY AND DISCUSSION 

In this paper we have considered the effect of the si- 
multaneous propagation of 2 spiral density waves moving 
at different pattern speeds, on the stellar radial velocity 
dispersion, an of a galactic disc. The numerical exper- 
iments performed concentrated on a primary 2-armed 
spiral wave and a weaker 4-armed one and their effect 
on the orbits of particles initially placed into circular 
orbits. This choice of spir al structure is motivated by 
recent COBE observations (Drimmcl & Spcrcel 2001). 

Following the growth of the spiral density waves, we 
find that the variance of the radial velocity resulting from 
2 spiral waves exceeds that expected from the sum of the 
variances from each pattern. In other words, the veloc- 
ity dispersion in an integration with both spiral density 
waves is larger than the sum of the effects of the individ- 
ual spiral density waves. The increase in <tr is particu- 
larly prominent at and near a corotation resonance with 
a primary, 2-armed spiral wave, in a combination with a 
secondary, 4-armed one at and near the OLR and ILR. 
Furthermore, we find that gr continues to increase near 
these resonances even after the spiral density waves have 
ceased to grow. This suggests that the coupling between 
the 2 spiral waves causes particles to diffuse in real space, 
even though the spiral density wave amplitudes are held 



fixed. 

We have examined the distribution of stars perturbed 
by one and two spiral density waves. When only one 
spiral density wave is present, little diffusion in phase 
space occurs. A comparison between these distributions 
(see Figures I9I13|) shows that diffusion is observed when 
multiple waves are present. This comparison confirms 
our hypothesis that the combined effect of the two spiral 
density waves causes stochastic diffusion of stars in phase 
space. 

Previously studied heating mechanisms have concen- 
trate d on heating from transient spiral density waves 
(ejr ISellwood fc Carlh^lT984l ICarlherg fc Sellwoodl 
19851 iDe Simone et alJ 120041 Uenkins fc Binnevl 119901: 
Sellw ood &: Binnevll2002|) . However in that case changes 
in the stellar velocity dispersion only occur during the 
growth or decay of spiral density waves. Here we observe 
an increase in the velocity dispersion even when the spi- 
ral density waves have ceased to grow. Hence, this is a 
different heating mechanism compared to those explored 
by these previous investigations. The heating we see here 
is similar to that observed at resonanc es in barred galax- 
ies. For example. iCombes et all ()1990j) observed stochas- 
tic heating in 3D simulations of barred peanut galaxies. 
We suggest that stochastic heating could also occur in 
disks of galaxies if more than one spiral density wave 
is present. In this case we could consider the source of 
the stochastic beh avior the overlap of the two resonant 
pertu rbations (e.g.. lMurrav fc Holmanlll999T : iLecar et alJ 
2001). 

We have established the possible combinations of pat- 
tern speeds for which heating is significant. This heating 
mechanism is strongly dependent on the relative veloc- 
ity between stars and the spiral pattern and thus on the 
galactocentric distance, R (e.g., see Figures 3,4). Conse- 
quently we expect strong variations in the velocity dis- 
persion as a function of position in the Galaxy. 

Future work could consider the role of Giant Molecu- 
lar Clouds (GMC) in addition to the currently explored 
mechanism by carrying out simulations in 3 dimensions. 
Here we have considered spiral density waves with fixed 
amplitudes and pitch angles. However future explo- 
rations could consider spir al structure consistent with 
Swing amplification (e.g., iFuchsl 120011) . Here we have 
focused on stars initially on circular orbits, however fu- 
ture explorations could consider stars born in spiral arms. 
Whereas in this paper we used a test -particle approach as 
was done bv lDe Simone et "all l)2004[) . future explorations 
can aim to produce self-consistent models where the or- 
bits of the stars are consistent with the spiral density 
waves themselves. 
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Fig. 1. — Radial velocity dispersion as a function of time for the primary 2-armed spiral pattern only. Each panel shows a different value 
of the angular velocity Q s ,a- The gravitational potential amplitudes are e a = —0.006 and e b = 0. The 2:1 Inner Lindblad Resonance (ILR) 
and Outer Lindblad Resonance (OLR) occur at Q s ,a = 0.3Slo and Q s ,a = 1.7f2o respectively. Time is in units of rotational periods and 
the radial velocity dispersion an is in units of the circular velocity Vq. Note the low values of an around the corotation radius Q s ,a = 1-0. 
There is no increase in an once the wave is grown. 
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Fig. 2. — Radial velocity dispersion as a function of time for the secondary, 4-armed spiral pattern only. Each panel shows a different 
value of the angular velocity Q 3: f The 4:1 Inner Lindblad Resonance (ILR) and Outer Lindblad Resonance (OLR) occur at f2 s i,=0.35,1.65, 
respectively. Since we incremented Q s (, by 0.1 this figure does not show a panel exactly at OLR or ILR. Thus we labeled the LRs on each 
side of their occurrence. Different line styles show different values of the spiral wave gravitational potential perturbation: = —0.001 
(dotted), e b = -0.002 (dashed), and e b = -0.003 (solid). . 
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Fig. 3. — Contours show Actr = <TR, a b — \l °\R,a 2 + °\R,f> 2 f° r 3 values of the amplitude of the secondary spiral wave, e;, = 
—0.001,-0.002,-0.003, and a constant value of the the primary amplitude, = —0.006 in (n SjQ ,f2 s j,)-space. Spiral amplitudes were 
grown simultaneously in 4 periods and were kept fixed after that. Here ctr j£1 and <jr : i, are the radial velocity dispersions as a result of 
only the primary or only the secondary spiral patterns, respectively. These were time averaged over 40 galactic rotations, excluding the 
time during which wave amplitudes grew (4 rotation periods). Similarly, ctr^i is the time average of a simulation run with both patterns 
perturbing the test particles. Consequently, Aor gives the additional heating caused by the combined effect of the two pcrturbers as 
opposed to the sum of the heating that would have been produces if each spiral wave propagated alone. In all panels Actr was normalized 
to its maximum value in the three panels allowing a proper comparison among the panels. The pattern speed of the primary wave, Q s>a , 
varies between 0.1 and 2.0 and that of the secondary one, Q a (,, varies between 0.5 and 1.6 in units of the angular rotation rate, Clo = 1 
at Ro = 1. The dashed lines show the locations of corotation and Lindblad resonances due to each of the two spiral wave perturbations. 
Next to each dashed line we labeled the resonances by denoting the corotation resonance with "CR", and the outer and inner Lindblad 
resonances by "OLR" and "ILR". For a given set of pattern speeds, a line at 45° slope predicts the variation of Aitr with radius, R, where 
the intercept of the line is set by the relative angular frequency, AU S = Q s ,a — £l s ,b- Large values of Aor are primarily observed along the 
lines Q s ,a = 1.0 and = 1.0, which correspond to the CR of each pattern, with the exception of the intersection of these lines; Aan(R) 
is a strong function of radius. 



12 



2:1 ILR 



CR 



2:1 OLR 



' .5 -eJ=-0.006, : e b =-0.001 - 

j j 4:1 OLR 



c 1 .0 



■+- 



CR 



0.5 



1 .5 -e»-0.006, : e b =-0.002 



•+- 



4:1 ILR 



4:1 OLR 



G 1.0- 



■+- 



CR 



0.5 
' .5 i 



+ 



4:1 ILR 



-0.006,: e b =-0.003- 

j j 4:1 OLR 



a 1 .0 



-+- 



CR 



0.5 

0. 



0.5 



4:1 ILR 



1 .5 



2.0 



Fig. 4. — Contour plots of the heating produced by the propagation of 2 spiral waves. Identical setup as in Fig. |3]but contours show 
ab> instead of Actr. Aside from the 2:1 LRs of the primary, 2-armed spiral wave, maximum values of er^ a (, are found at and around 
the CR with the primary and 4:1 LRs of the secondary patterns. As in Fig. |J] at the CR with both spiral waves very low radial velocity 
dispersion is seen (here it is a minimum). For all panels OR t ab was normalized to the largest value in the figure. 
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Fig. 5. — Time evolution of the radial velocity dispersion, an a i,, as a result of the simultaneous propagation of a 2- and a 4-armed spiral 
patterns (solid line). Time is in units of galactic rotations. The 2-armed spiral structure has the same parameters as the panel in FigUl 
for which Q s ,a = 1.0 (CR), whereas the 4-armed spiral wave has different angular velocity in each panel. The time evolution of the radial 
velocity dispersion due to the latter perturbation alone, an was shown in Fig. [5] For comparison, in addition to an a b (solid line) we 

also plotted <J an. a 1 + an,b 2 (dashed line) (see Eq. I1UI . The values of the perturbation amplitudes are e a = —0.006 and e;, = —0.001. We 
note the following dynamical effects around the 4:1 ILR and OLR in the case of an a (,: (1) an a j, increases up to a factor of m 3 compared 
to the dispersion expected from the noninteracting spiral waves, an, a' 2 + &R.,b 2 \ (2) an , a b grows with time; (3) resonances have effect 
over a larger range of angular velocities, compared to ^Jan, a 2 + <7R,b 2 i an d thus over a larger range of radii. 




Fig. 6. — Same as Figl^Jbut with the secondary gravitational potential amplitude tt = —0.002. The radial velocity dispersion at the 
LRs saturates faster than that in Fig. |S]as |ej,| was increased from 0.001 to 0.002. 




Fig. 7. — Same as Figures |2] and |g] but with the secondary gravitational potential amplitude tt = —0.003. Note the further decrease in 
the saturation time of an,ab as compared to Figures 151 and Ifil 
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Fig. 8. — Sa me as Figure 171 but with the primary spiral wave placing test particles at the 2:1 OLR (Q s ,a = 1.7), instead of the CR. As 
in Figures 15171 both spirals were grown simultaneously in the initial 4 rotation periods; the perturbation amplitudes are e a = —0.006 and 
e b = -0.003. 
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Fig. 9. — A snapshot of the evolution of a ring of stars subjected to a 2-armed spiral perturbation grown in the initial 4 rotation periods. 
The time is t = 40 rotations. The pattern speed is the same for each panel, f! s , a = 1-0, but the radius of the initial real-space particle 
distribution varies as: (a) R = 0.9, which places stars just inside the CR; (b) R = 1.0, placing stars exactly at the CR; and (c) R = 1.1, 
placing stars just outside the CR. The initial ring in each panel is indicated by the dashed line. The parameter a a has the default value of 
—6 which corresponds to a pitch angle p a ~ —18°. The curves on which stars lie formed at time t fx 4 rotation periods and remained fixed 
in the reference frame of the spiral structure for the rest of the run. 
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Fig. 10. — A snapshot of the evolution of a ring of stars subjected to a 4-armed spiral wave at and near the 4:1 ILR. The time is t = 40 
rotations as in Fig. [9] The pattern speed is the same for each panel, O s j = 0.65, but the radius of the initial real-space particle distribution 
varies as: (a) R = 0.9, which places stars just inside the 4:1 ILR; (b) R = 1.0, placing stars exactly at the 4:1 ILR; and (c) R = 1.1, placing 
stars just outside the 4:1 ILR. The initial ring in each panel is indicated by the dashed line. The parameter ct;, = —12, corresponding to a 
pitch angle pf, ~ —18°. The curves on which stars lie for R = 0.9, 1.1 were approximately constant in shape and orientation with time, as 
viewed in the rotation frame with the spiral wave. In contrast, exactly at the 4:1 ILR the rectangular curve spins clockwise with a period 
of about 9 rotations. 
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Fig. 11. — Same as Fig. 1101 but with the angular velocity of the secondary patt ern f2 3 f, = 1.35, which places stars at the 4:1 OLR at 
R = Rq. The pattern speed is the same for each panel, f! s ,a = 0.65. As in Fig. 1101 only when the initial particle distribution placed stars 
exactly at the resonance, panel (b), did the curve on which stars lie change shape and orientation with time. 
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Fig. 12. — Snapshots of the 3-step time evolution of 3 rings of stars subjected to both a 2-armed spiral wave near the CR and a 4-armcd 
one near the 4:1 ILR. This corresponds to the pa nel with Q s b = 0.7 in Fig. Time develops in each row from left to right. Note the very 
different morphology compared to Figures 1911 II Stars are observed to diffuse radially. All parameters of the primary spiral wave in each 
row are the same as those used in each row in Fig. [5] The difference here is that a second spiral wave is propagating in addition to the 
primary one and that we plot 3 steps of the time evolution since changes are taking place unlike in the time evolution of the systems shown 
in Fig. [5] Rows from top to bottom show particles subjected to 2 spiral wave perturbations moving at pattern speeds of Q s ,a = 1-0 and 
Q a b — 0.7 and initially positioned in a ring of radius: (a) R = 0.9, (b) R = 1.0, and (c) R = 1.1. Snapshots are at times i = 10,25,40. 
The dashed circles show where stars were positioned initially with a uniform azimuthal distribution. The crosses show the position of the 
galactic center. 
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Fig. 13. — Same as Fig. 1121 but with the secondary pattern speed near the 4:1 OL R, O s . t, = 1.2. This corresponds to the panel with 
f2 3 j, = 1.2 in Fig. [S] The radial diffusion of stars is notably increased compared to Fig. 1121 
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TABLE 1 
Parameters describing simulations 



Figures 


e a 


e 6 




° 6 


m a 


mi, 


n a 










-0.006 


-0.001,-0.002,-0.003 


-6 


-12 


2 


4 


0.1,0.2,. ..,1.9,2.0 


0.5,0.6,. ..,1.5,1.6 








-0.006 


-0.001,-0.002,-0.003 


-6 


-12 


2 


4 


0.1,0. 2,...,1. 9,2.0 


0.5,0.6,. ..,1.5,1.6 








-0.006 


-0.001 


-6 


-12 


2 


4 


1.0 


0.5,0.6,.. .,1.5,1.6 








-0.006 


-0.002 


-6 


-12 


2 


4 


1.0 


0.5,0.6,. ..,1.5,1.6 








-0.006 


-0.003 


-6 


-12 


2 


4 


1.0 


0.5,0.6,...,1. 5,1.6 








-0.006 


-0.003 


-6 


-12 


2 


4 


1.7 


0.5,0.6,. ..,1.5,1.6 








-0.006 


-0.002 


-6 


-12 


2 


4 


1.0 


0.65 


n 






-0.002 




-12 




4 




1.35 


12 




-0.006 


-0.002 


-6 


-12 


2 


4 


1.0 


0.7 






-0.006 


-0.002 


-6 


-12 


2 


4 


1.0 


1.2 



Note. - Spiral pattern parameters corresponding to simulations shown in the Figures. The 
perturbation strengths ei and ei are given in units of V 2 , the velocity of a star in a circular orbit at 
Rq. The pattern speeds, Q s ,a and Q s j,, are in units of Qq = Vo/Rq. The parameters a a and set 
the pitch angles of the spiral arms as m a cot(p a ) = a a and identically for m^jpi,, 05, where p a and 
Pi, are the pitch angles of the 2-armed and 4-armed spirals, respectively. Two different initial spatial 
distributions were used: (1) a uniform distribution in both R and tf> in the annulus (Ro — AR, Ro + AR) 
with An = 0.3 (Figures 11181 : (2) a uniform distri bution in the azimuthal direction and a delta peak 
at R = Rq, i.e., a ring of stars at R = Rq f Figures 191121 



